Anomalous conductivity: impact of nonlinearity and disorder 
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We reveal the intricate impact of nonlinearity and disorder on the thermal conductivity of acoustic 
chains. Disorder induces mobility edges and allows to control the amount of extended modes which 
are the ballistic channels for energy transport. Nonlinearity adds a diffusive conductivity channel 
through interacting localized modes, and controls the contact resistance at the edges. Analytical 
arguments and numerical results yield several crossovers between dominating heat channels, when 
varying the system size and the temperature. We demonstrate that the nonlinearity-induced inter- 
action between the modes alters the transport through existing conductivity channels and creates 
new ones, underpinning the observed phenomena. 
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Deviation from the normal (Fourier-law) heat conduc- 
tivity has recently become a matter of experimental and 
applied interest in the context of nanotube systems, in- 
spiring the research on solid-state thermal rectifiers and 
nanotube phonon waveguides The anomalous con- 
ductivity was clearly manifested in the pioneering exper- 
iments with carbon and boron-nitrid nanotubes, yielding 
K, oc N a , a~0.5-0.6[i|, with the other promise com- 
ing from nanowires In both systems phonons were 
reported to be the main heat flux carriers. Phonon scat- 
tering and interaction due to disorder (impurities) and 
nonlinearity (atomic interactions) are, thus, likely to de- 
termine the anomaly and its temperature dependence. 
Simulations with realistic models indicate that disorder 
Q or nonlinearity [U can give anomalous heat conduc- 
tivity, but the results on the exponent a are conflicting. 

One-dimensional momentum-conserving arrays serve 
as simple qualitative models to study anomalous heat 
conduction [||. Despite intense research and some quali- 
tative understanding, we still lack quantitative agreement 
on the main characteristics of anomalous conductivity 0- 
9J. Moreover, mostly harmonic chains with disorder, or 
anharmonic ordered chains were studied. Harmonic sys- 
tems do not equilibrate and the conductivity depends on 
the boundary conditions and the spectrum of thermal 
noise. For anharmonic ordered systems one lacks con- 
trol over the number of relevant long wavelength modes 
which contribute to anomalous conductivity. The inter- 
play between disorder and anharmonicity was touched in 
Refs. fiot 11 1, where the regime of normal conductivity 



[l(| was questioned due to finite size effects (ll| . 

In this Letter we uncover and study the intricate im- 
pact of nonlinearity on the thermal conductivity of the 
disordered Fermi-Pasta-Ulam (FPU) chain with fixed 
boundaries. Upon variation of the temperature and the 
chain size we observe transitions between the following 
regimes: (i) insulating behavior with a < 0, (ii) normal- 
like conductivity, a ~ 0, (iii) disorder-dominated anoma- 
lous conductivity, a ~ 0.52 — 0.58, (iv) nonlinearity- 
dominated anomalous conductivity, a ~ 0.38. The 



crossovers and the scaling of K. are explained by analyzing 
the properties and interaction of localized and delocalized 
modes, which are continued into the nonlinear regime, 
and emerging new heat conductivity channels. We show 
that anharmonic disordered systems offer a better way to 
study the mechanisms of anomalous conductivity, being 
also more realistic models for experimental setups. 

We consider the FPU-/3 chain of N equal masses, with 
disorder in the harmonic spring constants, and additional 
quartic anharmonicities in the spring potential: 



JV 



N+l 



4 ^n — 1 ) 



-(1 + Dn n )(x n - £„_i) 2 



(1) 



where x n (t) is the displacement of the n-th particle from 
equilibrium, p n (t) its momentum, K n G [—1/2,1/2] are 
random, uniform, and uncorrelated, (K n K m ) = a^S n ^ m . 
We apply fixed boundary conditions xq = xn+i = 0. 

The heat conductivity is measured by a standard ap- 
proach, when the thermal baths attached to the ends 
generate a temperature gradient and heat current along 
the chain 0|. We implement Nose-Hoover thermostats 
adding the terms — C±ii,JV to the respective equations of 
motion, where £± = x\ N /T± — l. The heat flux along the 
chain is defined as the time average of j = — | XX^n+i + 

n 

x n ){(l+DK n+1 )(x n+1 ~x n )+(3(x n+1 -x n ) 3 } pj. The heat 
conductivity coefficient reads K, = jN/ (T + — T_ ) then. 
We make use of the mean temperature T = (T+ + T_ )/2 
as a parameter corresponding to the energy density 
(E n ) = ksT = T, setting k B = 1 and (7+ - T_)/T = 0.5 
further on. 

We start with the canonical transformation to and 
from the normal modes of the harmonic lattice in the ab- 
sence of disorder: x n (t) — X) Q =i Qq(t) z qn which defines 
the mode space with N coordinates Q q (t) and eigenvec- 
tors z qn = ^/ jv+i sm ( N+\ ) • ^ n t ne presence of anhar- 
monicity and disorder the dynamics of the modes follows 
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Here uj q — 2 sin 2 (n+i) are ^ ne norma l mode frequencies. 
The coupling coefficients C q ^ p ^ yS [12| control the selective 
anharmonic interaction between modes, and the coef- 
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[13| reflect the all-to-all linear interaction due to dis- 
order. The nonlinearity and disorder parameters v — 
P/{N + 1), d = D/y/N+1 are small: v,d -C 1 for 
/3 = D = 1 and large system size TV ^> 1. 

Let us consider the case d ^ and z/ = 0. We aim 
at computing the new eigenvectors z qn defined through 
the transformation x n (t) = ^2 q —iQ q [t)z qn . We apply a 
perturbational approach for the harmonic mode go us- 
ing the small disorder parameter d: Q q (t) = Q q °\t) + 

dQq 1] {t) + . .., where Qq°\t) = for q ^ q . In the 
first order approximation, from Eq. ([2|) we obtain the 
equation of a forced oscillator for modes with q ^= q$: 

for the amplitude A of each mode 
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with the same scaling q c cx N 1 / 2 . These results tell, that 
for the disordered harmonic chain there is a thin layer of 
VTv" metallic modes in frequency space with frequencies 
< uj < uj c ~ iV -1 / 2 . For most practical purposes these 
modes can be described using the eigenvectors of the 
ordered harmonic chain. 

In a next step we add anharmonic terms, and compute 
the nonlinear analogs of normal modes - q-breathers [12] 
for the metallic modes of the disordered harmonic chain. 

We introduce {Q q ,P q } such that x n (t) — Q q (t)z qn . 

9=1 

It follows: 
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We develop the perturbation theory in powers of v, 
taking z qon as the zero order of approximation. Ne- 
glecting the higher order terms 0(isd), we obtain the q- 
breather [l2| with an energy distribution exponentially 
localized in the metallic mode space: 
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-u q u qo K q<qo Q q f . As a result we find About the mobility edge the modes get extremely weak 
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Accordingly, the time-averaged g-th mode energy is 



energies E qc ~ X qc ' qo E qo once A < 1 and the excitation 
of insulating modes is negligible. 

More than that, one can construct s-dimensional q-tori 
by continuing the first s -C q c disordered modes into the 
nonlinear regime, as done in the ordered case [l4| . These 
objects will be exponentially localized in q-space: 



{Eq) ~ 2 («s -J qo y ~ 2K - 4) 2 ■ (4) 

For Eq. ([3|) being valid we request (E qo +i) -C E qo , hence, 

go « g c = 2V2(N + l) 1/2 /Da K . (5) 

It follows that normal modes with mode 
numbers go -C q c approximately keep 
their plane wave eigenvector profile z qon = 



Ep/Eo cx (A T ) 2p , A T = pNE /s 



(8) 



where Eq and E p are the average energies of the modes 
1, s and (2p — l)s + 1, (2p + l)s [14]. Delocalization of 
tori happens when At > 1- Taking s cx q c we obtain the 
following strong stochasticity threshold: 
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real space in the presence of disorder. Therefore these 
metallic modes are still delocalized in real space. On 
the other side, Anderson localization implies that the 
eigenmodes of a one-dimensional disordered chain are 
localized. Therefore q c sets a mobility edge: for go <C g c 
the metallic eigenmodes are delocalized in real space, 
and for go > g c the insulating eigenmodes are localized. 
Note, that in short chains or/and at weak disorder 
all modes are metallic and delocalized if q c > N, ([5|) 
yielding N < 8/D 2 a 2 . Remarkably, the transfer matrix 
approach yields a lower boundary q c for spatially local- 
ized eigenstates (computing their localization length) Q 



Below the strong stochasticity threshold one can think 
of three basic heat conductivity channels (FigO}. (I) 
Metallic modes carry the heat flux ballistically by direct 
interaction with the heat baths, as in the linear case. 
The magnitude of interaction and the mode-specific heat 
flux are proportional to the square amplitudes next to 
the boundaries j q cx z 2 n=1N (boundary resistance) Q. 

For fixed ends j = ^2 q<Qc jq cx N~ 3 ^ 2 and K cx N~ x / 2 . 
(II) Insulating localized Anderson modes start to inter- 
act with each other due to anharmonicity, and are ex- 
pected to contribute a normal diffusive energy transport 
with /C(T) independent of N [l5(. (Ill) Anharmonicity 
also opens the third conductivity channel, by decreas- 
ing the boundary resistance for metallic modes. This 
happens due to thermalization among insulating modes, 
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FIG. 1: (Color online) Schematic representation of conduc- 
tance channels. Thick black metallic modes and thin black 
insulating modes are shown. (B)lue and (g)reen arrows show 
the energy fed into localized and delocalized modes channels 
by direct interaction with the heat baths. Both carry the heat 
fluxes (the latter more efficiently). Nonlinearity induces the 
heat flux via the localized modes into delocalized ones ((r)ed 
arrows), the third channel. 

which in turn will interact with the metallic modes in a 
broad range of space, thus eliminating the effects of fixed 
boundary conditions. 

The first conductivity channel (I) can be dominant 
only when the mean temperature T is low and the sys- 
tem length A is small enough. The heat current via the 
second localized mode channel (II) can be estimated us- 
ing recent studies of wave packet evolutions in disordered 
Klein-Gordon chains, where the harmonic approximation 
yields an Anderson insulator with a finite upper bound 
on the localization length of each eigenmode jlijj . For 
a given energy density (respectively finite temperature) 
the diffusion rate (and therefore also the heat conductiv- 
ity) are predicted to vary as K ~ T 2 for T > T cr and 
JC ~ T 4 for T < T cr . Here the crossover temperature 
T cr is a model-dependent temperature which depends on 
the strength of disorder. Notably this heat conductivity 
does not depend on the size of the system. Therefore 
the second channel will dominate over the first one if 
N > Ni2(T) where N%2(T) is some function which in- 
creases with decreasing temperature. The crossover to 
normal conductivity will be observed at A sa N\2(T). 

Similarly, the heat flux through the third conductivity 
channel can be estimated. Assume that 0(N) strongly 
localized insulating modes £ 9 = 0(1 ) are distributed uni- 
formly along the chain, E q — T, and heat the delocalized 
metallic mode go with initial energy E qo =0. It fol- 
lows Q qo +<Q 90 « /?£ P > 9c V P T 3 / 2 C(t) with C(t) be- 
ing an uncorrelated Gaussian noise with unit variance. 
Then (Q 2 ) cx f3 2 u; 2 T 3 t and the heat flux into the de- 
localized mode is J qo cx f3 2 u! 2 o T 3 . Extrapolating this 
result to the equilibrium conductivity problem, one ob- 
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FIG. 2: (Color online) Averaged heat conductivity coefficient 
vs. chain length in the log-log scale at different mean temper- 
atures and P = D = 1. The values of slope estimated by the 
linear fit are printed. Error bars estimate statistical error. 

tains J = . J 0n cx /3 2 T 3 ./V~ 1//2 and the associated 

heat conductivity K cx A 1 / 2 . Thus, for A > A 23 (T) we 
expect to observe a crossover from normal conductivity 
to anomalous a rs 0.5. Finally, with further increasing of 
the system size A, the strong stochasticity threshold ([9]), 
Eq = T, is reached for A > N sst . Then metallic modes 
start to strongly interact with each other, and theoretical 
predictions of fully developed turbulence within renor- 
malization group approaches, and mode coupling theo- 
ries may become applicable with the consequence that 
one will observe a ~ 1/3 .. . 2/5. 

To summarize, we expect that depending on the sys- 
tems size, heat is carried by: (i) ballistic metallic modes 
directly coupled to the heat reservoirs for A < A12 with 
a rs —1/2; (ii) insulating localized modes interacting 
with each other for A12 < A < A23 with a s» 0; (iii) bal- 
listic metallic modes which are coupled to the heat reser- 
voir via insulating localized modes for A23 < A < A sst 
with a « 1/2; (iv) strongly interacting metallic modes 
which are coupled to the heat reservoir via insulating lo- 
calized modes for N sst < A with a m 1/3... 2/5 . 

Let us turn to the numerical results, obtained with the 
integration time varying from 10 6 for the shortest chains 
to 2 • 10 7 for the longest, and up to 100 disorder real- 
izations. In all simulations /3 = D = 1. As observing 
all crossovers at a single temperature is currently im- 
possible (the required chain sizes and integration times 
are unattainable), we gradually increase T. The results 
need careful interpretation, as the finite size and residual 
temporal transient effects (even after long equilibration 
times) cannot be determined at present. 

At T = 10~ 4 we observe the first conductivity channel I 
with a < (Fig(2J). For larger temperature T — 0.001 we 
observe the normal-like conductivity of regime II in the 
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FIG. 3: (Color online) Size dependence of the heat conduc- 
tivity without and with additional random on-site harmonic 
potential (open and filled markers respectively). 

whole range of studied N, indicating that the localized 
modes conductivity channel is dominant. At T — 0.004 
we observe the crossover to the anomalous conductivity 
channel III with the value a ps 0.52 for large N, with even 
slightly bigger a at higher temperatures. It corresponds 
to the ballistic transport by metallic modes coupled to 
the heat baths via insulating localized modes. 

To confirm that this anomaly is due to the metallic 
ballistic modes, we add a random on-site harmonic po- 
tential 7(J2n=iQ ; n x ni fi ™ e [°! w o] to Q. Then a gap 
around zero frequencies opens in the frequency spectrum 
of the harmonic system, total momentum is not anymore 
conserved, and all modes become insulating and local- 
ized. Indeed, numerical results demonstrate the absence 
of anomalous behavior (FigJHJ) ■ The normal conductivity 
channel is affected too (due to increased localization of 
modes), but is still observable and becomes dominating. 

Finally, we observe the crossover due to the strong 
stochasticity threshold which yields a ps 0.38 (T = 0.2, 
FigJSJ, in which the temperature dependence of K be- 
comes indistinguishable. As all other crossovers, it can 
be observed for smaller values of TV with increasing T and 
at T = 4 this regime dominates starting from the small- 
est system sizes studied. Hence, the crossover N C (T) can 
be estimated as the intersection point of the linear fit for 
K(N) at T = 4 and the linear fit of JC(N) at a lower 
temperature, in its part where a ~ 0.5 is observed. For 
T = 0.02 and T = 0.2 one gets A c (0.02) ps 5.8 • 10 5 and 
N c (0.2) fa 4.5 • 10 3 , which agrees with N c cx T~ 2 © well. 

In conclusion, we shed light on the complex interplay 
between disorder and nonlinearity that determines the 
heat conductivity of acoustic chains, revealing a much 
more intricate picture than suggested before. Depending 
on the system size and temperature, one can find differ- 



ent conductivity channels to be dominant and observe 
crossovers as the parameters are changed. The mecha- 
nisms are: (i) ballistic transfer by metallic delocalized 
modes coupled directly to the heat baths, (ii) diffusive 
transfer by the insulating localized modes, (iii) ballistic 
transfer by metallic delocalized modes, the heat flux com- 
ing indirectly from the heat baths mediated via the insu- 
lating localized modes, (iv) turbulent transfer by metallic 
delocalized modes above the stochasticity threshold. The 
corresponding size dependence of the conductivity coeffi- 
cient is drastically different: insulating, normal, and two 
kinds of anomalous divergence. The studied system size 
is comparable to the number of atoms along nanotubes; 
therefore, the predicted crossovers may prove to be ob- 
servable even in the current experimental systems, if the 
temperature is varied. 
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